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ABSTRACT 


We  have  developed  a  random  walk  model  of  interstitial  diffusion  of 
light  impurity  atoms  in  a  host  lattice  of  heavy  atoms  in  the  presence  of  a 
thermal  gradient.  To  take  account  of  the  effect  of  the  thermal  gradient  on 
the  flux  of  impurity  atoms  we  introduce  a  bias  in  the  jump  direction  of  the 
interstitial  impurity.  We  assume  that  this  bias  is  due  to  the  temperature 
dependence  of  the  excluded  volume  effects  which  arise  during  the  jump  of  the 
impurity  atom  between  interstitial  sites.  The  resulting  random  walk  equation 
for  thermal  diffusion  is  consistent  with  both  positive  and  negative  heats  of 
transport  in  agreement  with  experimental  data.  Using  a  cell  model  approach 
and  the  assumption  of  local  equilibrium,  we  then  develop  equations  which  per¬ 
mit  the  calculation  of  the  bias  in  jump  direction.  The  theory  of  interstitial 
diffusion  developed  here  clarifies  and  supports  the  classic  Wirtz  model  for 
interstitial  thermal  diffusion. 


I .  INTRODUCTION 


The  dynamics  of  diffusion  in  a  crystalline  solid  maintained  at  con¬ 
stant  temperature  is  fairly  well  understood  at  the  present  time.  Both  the 

(1-3) 

simple  random  walk  models^  J  and  the  recently  developed  "dynamical" 

theories^  ^  give  much  insight  into  the  process  of  diffusion  and  produce  at 

least  qualitative  agreement  with  experiment.  For  the  more  complicated  problem 

of  diffusion  in  solid  subjected  to  a  temperature  gradient,  the  theoretical 

(7) 

description  is  less  satisfactory. 

(3) 

Allnatt  and  Rice  J  have  derived  general  random  walk  equations  which 
describe  the  flux  of  matter  in  a  crystalline  solid  with  both  concentration 
and  temperature  gradients  present.  To  apply  this  general  theory  to  any 
specific  problem,  certain  parameters  which  predict  both  the  magnitude  and 
the  direction  of  the  matter  flux  must  be  determined  outside  this  theory.  An 
example  of  such  a  parameter  is  the  heat  of  transport.  Using  the  random  walk 
model  together  with  certain  specific  physical  assumptions,  Allnatt  and  Rice 
concluded  that  the  heat  of  transport  is  always  positive  and  equal  to  the 
activation  energy  of  isothermal  diffusion.  Experiments,  however,  have  shown 
that  the  heat  of  transport  can  be  positive  or  negative,  depending  upon  the 
specific  system  under  study.  The  random  walk  model  itself  is  consistent  with, 
and  can  predict,  both  negative  and  positive  heats  of  transport.  Further  work 
on  the  determination  of  the  parameters  entering  into  the  random  walk  theory 
is  thus  desirable. 

In  order  to  be  able  to  examine  in  some  detail  the  physical  basis  of 
the  parameters  in  a  random  walk  model  for  thermal  diffusion  we  decided  to 
study  a  particularly  simple  system.  The  system  we  chose  is  that  of  inter¬ 
stitial  diffusion  under  the  influence  of  a  thermal  gradient  of  a  dilute 
solution  of  light  impurity  atoms  in  a  stationary  host  lattice  of  heavy  atoms. 


In  order  to  develop  the  relevant  random  walk  equations  we  consider  first 
interstitial  diffusion  in  an  isothermal  system  and  then  extend  our  equations 
to  interstitial  diffusion  in  a  thermal  gradient  (Sections  II  and  III).  By 
introducing  a  bias  in  the  jump  direction  of  the  random  walker  due  to  the 
thermal  gradient  we  are  able  to  obtain  both  positive  and  negative  heats  of 
transport.  To  evaluate  this  bias  we  introduce  a  simple  physical  model 
according  to  which  the  bias  is  due  solely  to  the  temperature  dependence  of 
excluded  volume  effects  which  arise  from  the  repulsive  forces  between  host 
lattice  and  impurity  atoms  during  the  jump  of  the  impurity  atom  between  inter 
stitial  sites  (Section  IV) .  This  simple  model  does  not  take  into  account 
other  possible  mechanisms  which  might  also  produce  a  bias  in  jump  direction 
under  the  influence  of  a  thermal  gradient  such  as,  for  instance,  chemical 
host-impurity  interactions.  From  our  model  we  then  calculate  the  jump  bias 
to  a  first  approximation  and  use  it  to  estimate  heats  of  transports  for  the 
diffusion  of  H  and  C  in  a-(bcc)Fe.  These  estimated  heats  of  transport  are 
in  order  of  magnitude  agreement  with  experimentally  determined  ones  and  have 
the  correct  sign. 

We  find  that  there  is  a  close  relationship  between  the  random  walk 
approach,  Rice's  "dynamical"  theory ^  and  the  Wirtz  model ^  for  thermal 
diffusion.  These  relationships  are  discussed  in  some  detail  in  the  last 
section  of  this  paper. 


II.  INTERSTITIAL  DIFFUSION  OF  IMPURITIES 
A.  Isothermal  Diffusion 

In  order  to  describe  the  isothermal  diffusion  of  impurity  atoms  via 
an  interstitial  mechanism  we  make  the  following  assumptions:  (a)  the  light 
impurity  atoms  diffuse  in  a  stationary  host  lattice  of  heavy  atoms,  (b)  the 
impurity  atom  concentration  is  sufficiently  low  that  interstitial  lattice  sites 
on  either  side  of  a  given  interstitial  impurity  atoms  can  be  taken  as  unoc¬ 
cupied  and  that  we  can  work  with  concentration  rather  than  activities,  i.e. 
the  system  is  ideal,  (c)  the  diffusion  is  one  dimensional  (in  cubic  crystals 
diffusion  in  the  three  normal  directions  can  be  treated  independently). 

Consider  now  two  interstitial  planes  1  and  3  which  are  a  distance  a 
apart.  (See  Figure  1.)  Let  the  number  of  impurity  atoms  per  unit  area  on 
plane  1  be  n(l)  and  on  plane  3  be  n(3).  The  total  jump  rate  I\  which  is 
assumed  to  be  a  function  of  temperature  only,  is  defined  so  that  n(l)T(l)  gives 
the  total  number  of  atoms  per  second  per  unit  area  jumping  from  plane  1  to 
planes  on  both  sides  of  1.  We  wish  to  consider  the  net  flux  J  of  impurity 
atoms  across  the  host  lattice  plane  2  between  interstitial  planes  1  and  3. 

This  is  clearly  given  by 

J(2)  =  PR(l)ni)n(l)  -  PL(3)I\3)n(3)  (2J) 

where  PR(1)  represents  the  probability  on  a  given  jump  that  an  atom  on  plane 
1  will  jump  to  the  right  and  where  P^(3)  is  the  probability  of  an  atom  on 
plane  3  jumping  to  the  left.  Right  and  left  refer  respectively  to  jumps 
toward  positive  or  negative  x  in  Figure  1.  In  the  isothermal  case  jumps  to 
the  right  and  left  are  equally  probable  so  PR(1)  =  P^(3)  =  J  and  t*ie  j^P  rate 
IXT)  is  the  same  in  planes  1  and  3.  Equation  (2.1)  then  becomes 

J(2)  =  ^jn(l)  -  n(3)  ]  IXT) 

-  4  - 


(2.2). 


ac(l),  where  c(l)  is  the  number  of  atoms  per  unit  volume  around 

1  nnn  1 


plane  1  we  can  write 


J(2)  =  j  a  [  c(l)  -  c(3)  J  rCT) 


'  1  (2.3) . 
We  assume  that  the  concentration  is  a  slowly  varying  function  '  x  over  the 
range  of  the  jump  distance  a  and  write 


c(l)  -  c(3)  =  -a 


Equation  (2.3,  now  takas  the  for.  of  Pick's  ,irst  la„  o£  diffusion; 


J  =  -  J  a2r(T)  |£ 


(2.4). 


(2.5) 


We  can  then  identify  the  emnirimi  • 

empirical  diffusion  constant  D(T)  in  Fick's  first 


J  =  -D(T) 


(2.6) 


D(T)  =  I  a2r(T) 


^  (27 

Ihus,  as  is  well  known,  the  isothermal  jump  rate  T(T)  can  be  determined  from 
measurements  of  the  diffusion  constant  D(T)  over  a  range  of  temperatures. 


~ — -lffusion  j£  a  Temperature  Gradient 

We  consider  now  the  effect  of  imposing  a  temperature  gradient  across 

C"Stal-  ThC  "°St  dif£usi0”  -  to  be  negligible  ana  . 

well-defined  texture  distribution  T(x,  is  therefore  established  which  can 

*  defl"e<‘  ln  °f  the  oscillations  of  the  host  lattice  atoms 

around  their  equilibrium  positions.  a  given  experimental  situation,  this 

temperature  gradient  could  either  be  measured  by  appropriate  probes  or  calcu- 

ated  by  Fourier's  law.  A  local  temperature  a  Tfxp  and  local  temperature 


5 


(2.8). 


gradient  AT(Xj)  at  plane  1  can  be  defined  by 


ATCxp  s  ATj 


&T(x) 

dx 


x=x. 


T(x3)  -  T(Xj) 


We  assume  that  the  total  jump  rate  r  for  the  impurity  atoms  on  plane 
1  in  the  presence  of  a  temperature  gradient  AT  is  a  function  of  and  AT^  such 
that  when  ATj=0  the  jump  rate  is  that  of  an  isothermal  system  at  temperature 
Tj .  We  expand  r  in  a  Taylor  series  in  AT 

IXT,AT)  =  T(T)  +  A (T) AT  +  B(T) (AT) 2  +  ...  (2.9) 


where  T(T)  is  the  total  jump  rate  in  an  isothermal  system  and  A(T)  and  B(T) 
are  functions  describing  the  effect  of  the  temperature  gradient  on  the  jump 
rate  I\  However,  the  term  proportional  to  AT  must  vanish,  i.e.  A(T)=0, 
because  the  total  jump  rate  T(T,AT)  must  not  change  if  the  direction  of  the 
temperature  gradient  is  reversed.  Thus  to  first  order  in  AT,  the  total  jump 
rate  I\T)  at  a  local  temperature  T^  and  a  local  temperature  gradient  AT^  is 
the  same  as  that  in  an  isothermal  system  at  the  temperature  T^.  The  tempera¬ 
ture  gradient,  to  first  order,  thus  does  not  cause  a  change  in  the  total 
jump  rate  but  rather  a  bias  in  the  jump  direction  --  i.e.  a  redistribution  of 
the  fraction  of  molecules  going  to  the  right  or  the  left.  For  the  isothermal 
case,  PR(T,AT=0)  =  P^(T,AT=0)  *  ■j;  in  a  temperature  gradient  we  write  an  ex¬ 
pansion  in  AT 

PR(T,AT)  =  j+  W  (T)  AT  +  ...  (2.10) 

where  W(T) AT  gives  the  bias  in  the  jump  direction  to  first  order  in  AT  due  to 
the  effects  of  the  temperature  gradient  on  the  host  lattice.  Analogously 

PL(T,AT)  =  I  -  W(T)AT  v  ...  (2.11), 

where  we  assume  the  temperature  increases  in  the  positive  x  direction.  From 
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Eqs .  (2.1),  (2.10)  and  (2.11)  wc  now  find  for  the  flux  across  plane  2 

J(2)  =  [y  +  WCT^AtJ  rCrpnd)  -  [y-  K(T3)AT3]  r(T3)n(3)  (2.12). 

Note  that  the  IXT)  refer  to  the  jump  rates  in  an  isothermal  system  (at  tempera 
tures  Tj  and  T3  respectively)  and  that,  to  first  order,  the  effect  of  the 
temperature  gradient  shows  up  only  in  the  W's.  We  take  ATj  =  AT3  since  the 
temperature  gradient  cannot  change  significantly  over  distances  of  the  order 
of  a  lattice  spacing.  Also  since  W(T)AT  is  already  first  order  in  AT  we  can 
set  W(T3)AT  =  W(Tj)AT  in  Eq.  (2.12)  to  obtain 

J (2)  =  [y  +  WCT^AtJ  rd^nd)  -  [y  -  WCT^ATj]  r(T3)n(3)  (2.13). 

Setting 

n (3)  =  n(l)  +  |y  a  (2.14) 

r(T3)  =  HTj)  +  ^|P-AT  (2.15), 

substituting  the  relations  (2.14)  and  (2.15)  into  Eq.  (2.13)  and  keeping  only 

terms  proportional  to  the  first  power  of  AT  or  ,  we  obtain 

i  ;wn  i  am\) 

J (2)  =  -  yarerp  ^IILL-  y— y~-ATn(l)  +  2W(Tj) ATnTj)n(l)  (2.16). 

The  use  of  the  relation  n^  =  ac^  and  Eq.  (2.8)  finally  yields  the  thermal 
diffusion  equation 

J (2)  =  -  y  a2r(T)  ||  -  y  a2c!tT)  |£  -  4W(T)J  (2.17) 

he  ftT 

where  nT>>  sr  •  c-  srand  w  are  all  to  be  evaluated  on  plane  1.  Using  Eq. 
(2.7),  which  relates  r(T)  to  the  isothermal  diffusion  constant  D(T),  we  can 
also  write  Eq.  (2.17)  in  the  form 

J (2)  =  -D(T)  |  |y  -  c  |I  . 4W(T)J  |  (2.18) 

The  first  term  in  Eq.  (2.18)  is  the  familiar  isothermal  result.  The  second 
term  which  is,  proportional  both  to  the  temperature  gradient  and  the  local 


concentration,  reflects  the  effect  of  the  temperature  gradient.  Our  sub¬ 
sequent  discussion  is  concerned  with  the  analysis  and  interpretation  of  the 


thermal  diffusion  F.qs.  (2.17)  and  (2.18). 

We  can  relate  our  jump  bias  parameter  W  to  an  experimental  parameter 
Q^,  usually  referred  to  as  the  "measured  heat  of  transport."  It  is  defined  by 
the  equation 


-kT 


dftnc 

dx 


n  dfcnT 
Sn  dx 


(2.19) 


from  which  is  to  be  evaluated  when  the  measured  concentration  gradient 
of  the  diffusing  species  has  become  independent  of  time.  In  our  approach  the 
concentration  gradient  becomes  independent  of  time  when  there  is  no  net  flux. 
Thus  we  set  J(2)  =0  in  Eq.  (2.18)  and  rewrite  the  resulting  equation  as 

dinT 


-•‘t2§T--  ^(tT  -  4W) 


(2.20). 


Comparing  Eqs.  (2.19)  and  (2.20)  we  sec  that  our  random  walk  theory  relates 
W  to  by 

r2  (  dfcnD 


kT 


(Ir  •  «) 


(2.21). 


It  will  be  noted  that  the  heat  of  transport,  is  given  by  the  difference 
between  two  quantities.  The  first  of  these,  can  be  evaluated  from 

measurements  on  isothermal  systems.  The  second,  the  jump  bias  W(T),  can  be 
estimated,  in  principle,  on  the  basis  of  an  appropriate  physical  model.  We 
present  such  a  model  in  Section  IV.  It  is  clear  from  Eq.  (2.21)  that  can 
be  positive  or  negative  depending  upon  the  relative  magnitudes  of  and  W. 
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111.  PHYSICAL  INTERPRETATION  OF  THE  THERMAL  DIFFUSION  EQUATION  (2.18) 


Equations  (2.17)  and  (2.18)  give  tlic  flux  across  a  lattice  plane  as  a 
function  of  the  diffusion  constant,  the  concentration  gradient  and  the  tempera¬ 
ture  gradient.  The  only  term  in  (2.18)  which  cannot  be  evaluated  from 
isothermal  measurements  is  W,  the  bias  in  jump  direction  due  to  the  effect  of 
the  temperature  gradient  on  the  host  lattice.  An  a  priori  evaluation  of  W  is 
very  difficult.  It  could  be  taken  as  a  disposable  parameter  to  be  evaluated 

from  experimental  data.  This  approach  would  be  similar  to  that  taken  by 
(9) 

Girifalco  and  our  theory  can  easily  be  related  to  his. 

We  shall  now  indicate  how  our  model  of  interstitial  thermal  diffusion 

and  the  flux  Eqs.  (2.17)  and  (2.18)  resulting  from  it  can  be  used  to  interpret 

the  available  experimental  data.  It  is  clear  that  or 

ol 

tivc  quantity  —  the  total  jump  rate  T  certainly  increases  as  the  temperature 
is  increased,  and  experimental  measurements  confirm  that  the  diffusion  con¬ 
stant  is  an  increasing  function  of  the  temperature.  We  would  also  expect,  in 
general,  that  W  is  positive  since  a  temperature  gradient  biases  the  jump 
direction  towards  the  region  of  higher  temperature  (which  for  simplicity  we 
shall  call  the  hot  side).  This  is  due  to  the  fact  that  the  host  lattice  atoms 
have,  on  the  average,  larger  vibration  amplitude?  on  the  hot  side  and  thus 


is  a  posi- 


thc  interstitial  atom  can  jump  more  easily  in  this  direction.  Equation  (2.20) 
shows  that  it  is  the  interplay  between  W  and  — that  determines  the  steady 
state  concentration  gradient  established  by  a  temperature  gradient.  If  the 
lattice  bias  effect  is  sufficiently  large  so  that  4W  >  SillS  then  ~  >  0  in  the 
steady  state,  i.c.  there  is  a  net  motion  of  the  impurity  atoms  towards  the  hot 


side  starting  from  an  initial  homogenous  distribution.  When  the  lattice  bias 

K  is  sm'.ill  and  4W  <  then  ^-<0.  In  this  case  there  is  a  net  motion  of 

oi  ox 

impurity  atoms  towards  the  cold  side.  This  would  tend  to  be  the  ease  for  small 
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interstitial  impurities.  In  the  limit  of  no  lattice  bias  on  the  jump  direction, 

dc 

i.c.  for  W=0,  one  lias  ^-<0  and  the  impurity  atoms  always  move  toward  the  cold 

side  simply  because  the  total  jump  rate,  ITO  is  higher  in  the  high  temperature 

region  than  in  the  low  temperature  one. 

(3) 

Allnatt  and  Rice  implicitly  assumed  K=0  in  their  random  walk  treat¬ 
ment  of  interstitial  diffusion  and  thus  concluded  that  is  always  positive 
and,  from  Eq.  (2.21),  equal  to  the  activation  energy  for  isothermal  diffusion, 
kl?  t  jn  the  gCncrai  casc  where  there  is  a  lattice  bias  effect  one  must, 

however,  consider  both  terms  in  Eq.  (2.21)  to  predict  the  sign  and  magnitude 
of  the  heat  of  transport. 
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IV.  ESTIMATION  OF  THE  LATTICE  BIAS  EFFECT 
A.  Isothermal  Cell  Model 

We  assume  that  the  lattice  bias  effect  W  is  primarily  determined  by 

the  relative  size  of  the  impurity  atom  as  compared  to  the  interstitial  cavity  , 

in  the  host  lattice.  The  impurity  atom  causes  a  lattice  strain  and  a  jump 

from  one  interstitial  site  to  another  can  occur  only  when  (a)  the  amplitude 

of  vibration  of  the  impurity  atom  is  large  and  directed  towards  the  adjacent 

interstitial  site,  and  (b)  the  atoms  immediately  surrounding  the  impurity 

atom  move  apart  to  permit  the  impurity  to  fit  between  them,  and  there  is  at 

the  same  time  a  local  expansion  of  the  lattice  at  the  adjacent  interstitial 

site  which  permits  the  impurity  to  fit  comfortably  into  its  new  site.  This 

(4) 

dynamical  description  of  diffusion  has  been  discussed  extensively  by  Rice. 

In  a  temperature  gradient  it  seems  reasonable  that  condition  (b)  occurs  more 
frequently  on  the  hot  side,  where  the  lattice  atoms  will  more  often  have 
sufficient  energy  to  reach  the  "critical  configuration."  This  then  establishes 
a  bias  in  the  jump  direction. 

We  present  below  a  method  for  the  calculation  of  W  which  makes  use  of 
the  ideas  in  Rice's  dynamical  theory  of  diffusion.  This  treatment  also  leads 
to  a  more  precise  description  of  the  ideas  behind  the  Wirtz  model'  '  and  sug¬ 
gests  that  it  has  some  validity  in  the  description  of  interstitial  diffusion. 

C4) 

Following  Rice,  J  we  assume  that  a  particle  moves  irreversibly  to 
a  neighboring  interstitial  site  whenever  (a)  fluctuations  cause  it  to  achieve 
a  critical  displacement  along  the  line  joining  the  old  to  the  new  site  and  (b) 
all  host  lattice  particles  which  might  interfere  with  the  jump  have  moved  out 
of  the  way.  The  jump  rate  r  is  given  by  the  frequency  with  which  conditions 
(a)  and  (b)  occur  simultaneously.  Rice  then  makes  the  physically  reasonable 
approximation  that  condition  (a),  which  involves  a  very  large  amplitude 
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fluctuation,  occurs  less  frequently  than  (b) .  The  jump  rate  in  a  particular 
direction  X  can  then  be  written  symbolically  as 

1^  =  P^r=  v(a)p(b)  (4.1) 

where  v(a)  is  the  frequency  with  which  condition  (a)  occurs  and  p(b)  is  the 
probability  that  the  neighboring  host  lattice  atoms  have  the  appropriate 
configuration  for  the  passage  of  the  interstitial  impurity  independent  of  the 
position  of  the  jumping  interstitial  atom.  In  our  one  dimensional  model  the 
total  jump  rate  T  is  given  by  +  1^* 

The  computation  of  v(a)  is  very  difficult  since  the  motion  of  the 
interstitial  impurity  during  a  jump  involves  large  fluctuations  which  certainly 
cannot  be  well  described  by  the  harmonic  approximation  customarily  used.  Further 
more  the  critical  amplitude  of  this  fluctuation  leading  to  diffusion  and  the 
appropriate  "reaction  coordinate"  cannot  be  obtained  from  the  dynamical  theory 
itself.  Fortunately,  under  certain  approximations  to  be  detailed  below,  the 
lattice  bias  effect  in  a  thermal  gradient  depends  only  upon  p(b)  so  that  one 
does  not  have  to  calculate  v(a). 

Rice  has  pointed  out  that  in  the  isothermal  case  the  probability  p(b) 
could  be  estimated  fairly  accurately  by  a  'harmonic  potential  model  since  the 
vibrations  of  the  host  lattice  atoms  will  be  of  relatively  small  amplitude 
about  their  equilibrium  positions  and  since  the  equilibrium  positions  of  the 
host  lattice  atoms  are  fixed  in  space.  To  extend  the  calculation  of  p(b)  to  a 
crystal  in  a  thermal  gradient  we  make  two  basic  assumptions.  First,  we  assume 
a  cell  model  for  the  crystal  and  second,  we  assume  that  the  magnitude  of  the 
vibrations  of  a  host  lattice  atom  in  its  cell  is  a  function  of  the  local  cell 
temperature.  In  addition  to  these  fundamental  approximations  we  will  make 
several  additional  non-essential  approximations  for  computational  and  concep¬ 
tual  simplicity.  We  will  first  discuss  the  calculation  of  p(b)  for  an 
isothermal  system. 
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Since  most  thermal  diffusion  experiments  occur  at  high  temperatures 
(500°C  to  2000°C) ,  it  seems  reasonable  to  consider  the  simple  Einstein  inde¬ 
pendent  oscillator  model  to  describe  the  host  lattice  vibrations .  To 

apply  this  model  to  our  calculation  of  p (b)  we  make  use  of  a  single  occupancy 
cell  model  for  the  crystal  (in  the  presence  of  the  impurity  atoms)  where  each 
atom  moves  within  its  cell  in  the  potential  fields  of  its  neighbors  which  are 
taken  to  be  fixed  at  their  static  equilibrium  positions.  This  is  simply  the 
Lennard-Jones  and  Devonshire  theory,  which  though  developed  for  liquids, 

fl  2) 

is  really  a  theory  for  the  solid  state.  To  simplify  the  problem,  the 

nearest  neighbors  are  now  treated  as  uniformly  "smeared”  over  a  spherical 
surface.  This  yields  a  spherically  symmetric  potential  which  can  be  expanded 
to  second  order  about  the  potential  minimum  and  in  this  form  leads  to  the 
Einstein  independent  oscillator  model.  Let  r^  represent  the  displacement  of 
molecule  i  from  its  equilibrium  position  at  the  center  of  its  cell.  The 
subscript  1  will  refer  to  the  interstitial  atom.  For  this  simple  cell  model 
the  classical  canonical  n  particle  probability  density  (r-y^. .  «£)|+]^  .  that 
is,  the  probability  density  for  finding  host  lattice  atom  2  in  dj^  at  atom 
3  in  d^  at  etc.,  is  the  product  of  single  probability  densities 

=  p(1)Cr2)P(1)(r3)...PCi;>(rn+1)  (4.2) 


where 


,(D 


-BV.(r.) 


fe  1 


(4.3) 


dr. 

~i 


and  P  =  pjr  .  Here  V^(r^)  is  the  spherically  symmetric  (smeared)  potential 
which  atom  i  feels  when  displaced  a  distance  r^  from  its  equilibrium  position 
Within  the  harmonic  lattice  approximation,  we  expand  V^(r^)  about  thi 
minimum  r.=0  to  obtain 

l 


V.  (r.)  =  V. (0)  +  i  K.r? 

iv  i'  iv  '  2  l  l 


(4.4). 
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We  furthermore  assume  that  the  normalization  integral  in  Eq.  (4.3),  which 
properly  should  extend  only  over  the  i™  cell,  can  be  extended  to  infinity 
since  V^(r^)  dies  off  so  rapidly.  With  these  approximations  Eq.  (4.3) 
becomes 


,d) 


<v  ■& 


«i  2 

3/2  -^r. 


(4.5) 


The  probability  density  (£2x,y  •  •JEn+p  in  Ecl-  (4-2)  can  now  be  calculated 
for  all  those  n  host  lattice  atoms  which  can  conceivably  interfere  with  the 
jump  of  the  interstitial  atom  in  a  particular  direction. 

The  physical  basis  of  our  calculation  is  that  a  jump  can  take  place 
only  when  the  host  lattice  atoms  have  moved  in  such  a  way  that  a  passage  has 
been  opened  for  the  interstitial  atom.  Note  that  by  considering  only  those 
particular  configurations  of  the  host  lattice  atoms  which  permit  a  jump  we 
are  in  effect  improving  the  inaccurate  treatment  of  the  repulsive  forces 
between  the  host  lattice  and  interstitial  atoms  inherent  in  the  harmonic 
lattice  model. 

We  assume  for  simplicity  that  the  repulsive  forces  between  the  host  and 
interstitial  atoms  can  be  approximated  by  hard  core  interactions.  The  host  lattice 
configurations  which  permit  a  jump  of  the  interstitial  atoms  can  then  be  estimated 
from  the  geometry  of  the  system.  From  the  geometry  one  can  determine  a  critical 
displacement  of  the  i*^  lattice  atoms  along  a  particular  direction  ^  which 

C41 

is  essential  for  the  passage  of  the  impurity  between  interstitial  sites.  1 
For  example,  the  nearest  neighbor  cell  atoms  to  the  impurity  must  move 
sufficiently  far  in  a  direction  perpendicular  to  the  "jump  path"  of  the 
impurity  that  the  radius  of  the  resulting  hole  is  equal  to  or  larger  than  the 
hard  sphere  radius  of  the  interstitial  atoms.  Motions  of  these  host  lattice 
atoms  in  other  directions  are  probably  less  important  and  for  simplicity,  we 
assume  that  they  do  not  effect  the  jump  probability.  With  these  assumptions, 
whose  validity  can  admittedly  be  verified  only  a_  postcrori,  the  probability 
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p(b)  can  be  written  as  the  product  of  the  probabilities  JL(£^)  for  all  the  n 
neighboring  atoms  that  the  1  *  host  lattice  atom  lias  a  displacement  along 
equal  to  or  greater  than  the  necessary  critical  displacement  From  Eq. 


(4.5)  we  find 

W  =  \  erfe 


(4.6) 


where  the  integration  over  r^  is  from  to  infinity  in  the  direction  )L  and 
from  minus  infinity  to  infinity  in  the  other  normal  directions.  The  comple¬ 


mentary  error  function  used  in  Eq.  (4.6)  is  defined  by 


erfc  z 


■  fe  f  •" 


(4.7). 


Kith  these  approximations,  the  jump  rate  TR  of  our  one -dimensional  model  in 
the  absence  of  a  temperature  gradient  then  becomes  [see  Eq.  (4.1)] 


rR(TrAT=0)  =  v(a)p(b)  =  vCa.Tj)  17  O.^.Tj) 


(4.8) 


where  we  have  explicitly  indicated  the  temperature  dependence  of  the  terms. 


AX^  is  the  distance  from  plane  1  at  temperature  to  cell  i  at  temperature 
and  £  is  the  jump  length  defined  earlier.  Thus,  in  a  temperature  gradient, 
Eq.  (4.8)  becomes 
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rR(TrAT)  =  vCa.Tj)  77^  fi.CC./fj+a.AT) 


(4.11). 


The  frequency  v(a)  in  our  local  equilibrium  model  is  a  function  only  of  the 
local  temperature  Tj  of  the  cell  around  the  impurity 'atom  and  not  of  the 
temperature  gradient.  It  is  now  easy  to  relate  the  quantities  of  Eq.  (4.11) 
to  W(T),  the  parameter  in  our  random  walk  model.  By  definition,  to  first 
order  in  AT,  we  have 

1  rR(T  AT) 

Pr(T1,AT)  =  y  +  W(Tj) AT  =  rR(Ti,AT)+rL(Ti,AT) 


rR(T,,AT)  f  d£nrR(T  ,AT=0)  1 

"  2  rR(T1,AT=0)  =  2  +  dTATl  ATJ  (4,12) 

where  we  have  used  the  fact  [see  Eq.  (2.9)]  that  the  total  jump  rate  rR  + 
is  the  same,  to  first  order  in  AT,  as  that  for  the  isothermal  system  with 


(4.12) 


AT=9.  Substitution  of  Eq.  (4.11)  into  (4.12)  yields  for  W(Tj) 

n+1  3£nn.U.,V) 

- 1  £ — iar-L-“i  '4-13>' 

Note  that  v(a,Tj),  which  would  be  difficult  to  calculate,  drops  out  of  the 
theory  as  discussed  above.  The  bias  in  jump  direction  is  thus  due  to  the 
different  probabilities  p(b)  with  which  the  host  lattice  can  achieve  the 
critical  configurations  on  either  side  of  the  impurity  in  a  thermal  gradient. 

Equation  (4.13)  allows  us  to  relate  the  random  walk  parameter  W  to 
the  local  lattice  geometry  and  the  force  constants  K. .  From  Eq.  (4.6)  we 


obtain 


hi  \  -n 

■jKj.y  W 


j; 


where 


(h<r 

"  S  \2^7/ 


(4.14) 


(4. IS) 
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, .  •  (13) 

Using  the  inequalities 


,  2  ZTfi 

z+(z  +2) 


r  <  c 


F 


dt  *  — /  2  4  \l'/7 

T  * ") . 


(4.16), 


Eq.  (4.14)  becomes 


9X,nf2^  (C  ^  >1  j)  y  j 


where 


dT 


1  ^  n  <  2 
u 


2T, 


[i *  o.n/r2),/2] 


(4.17) 


(4.18) 


Usually  y2  »  n  since  \  KjE2,  the  potential  energy  of  the  host  lattice 
oscillator  at  its  critical  displacement,  is  much  greater  than  Wj.  the  average 
energy  per  degree  of  freedom  of  the  host  lattice  oscillator  at  temperature  T,. 
To  a  first  approximation  we  can  therefore  neglect  the  n/y2  term  in  Eq.  (4.17) 

and  write 


(4.19). 


&T 


2kT 


We  then  finally  obtain  from  Eq.  (4.13),  for  the  bias  W(Tj)  in  the  jump  direc- 
tion 

1  »S.WKiCiV  (4.20). 

W(T1>  -  JTJ-  S2  ^  )  “i 

The  lattice  bias  around  the  impurity  site  'l  is  thus  simply  related  to  the  sum 
of  the  reduced  potential  energies  of  the  n  nearest  neighbor  independent 
oscillators  (atoms)  of  the  host  lattice  at  their  critical  displacement 

weighted  by  the  distance  parameters  ou. 

If  we  assume  that  the  isothermal  diffusion  constant  can  be  written 


as 


D(T)  »  Doe 


AE 

FT 


(4.21) 


where  AE  is  the  activation  energy  for  diffusion,  then  the  heat  of  transport, 
defined  in  Eq.  (2.21),  becomes,  using  Eqs .  (4.20)  and  (4.21), 
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(4.22). 


The  heat  of  transport  for  interstitial  thermal  diffusion  can  thus,  in  principle, 
be  calculated  from  Eq.  (4.22)  if  the  isothermal  activation  energy  for  diffusion 
is  known. 

There  are,  however,  some  ambiguities  in  calculating  the  parameters. 

JL,  and  for  any  given  host  lattice -impurity  system.  These  ambiguities 

are  compounded  by  the  paucity  of  experimental  data  for  Qm  for  the  type  of 

interstitial  diffusion  we  consider  here  —  reliable  data  seems  to  be  available 

(7) 

only  for  a  few  systems.  This  lack  of  data  and  the  ambiguities  inherent 
in  the  calculation  of  K^,  and  cu  makes  it  impossible  to  check  at  this  time, 
by  comparing  calculated  and  measured  heats  of  transport,  whether  heats  of 
transport  calculated  from  Eq.  (4.22)  are  in  reasonable  agreement  with  experi¬ 
mental  values.  We  have  obtained  order  of  magnitude  agreement,  as  well  as  the 

(7a) 

correct  sign,  for  for  the  thermal  diffusion  of  C  and  H  in  a-(bcc)Fev  J  by 

(14) 

calculating  from  compressibility  data  using  the  simple  Einstein  theory 
and  estimating  and  n  from  the  atomic  radii  of  C,  H  and  Fe  and  the 

known  crystal  structure  of  the  body  centered  cubic  Fe  lattice.  This  agreement 
is  encouraging  and  suggests  that  this  theory  should  get  a  more  careful  test 
when  additional  data  becomes  available.  Hopefully,  this  attempt  at  a  theory 
will  stimulate  some  relevant  experimental  work. 
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V.  RELATION  TO  THE  WIRTZ  MODEL 


The  spirit  of  this  treatment  and  the  resulting  liq.  (4.22)  are  very 
reminiscent  of  the  classic  Wirtz  model. ^  The  Wirtz  model  for  interstitial 
diffusion  can  be  interpreted  to  state  that  the  jump  rate  to  the  right,  1^,  is  . 
given  by 

rR  =  vexp  irfwWV1*)  ^A) 

where  v  is  a  vibration  frequency.  The  quantity  E^  is  interpreted  as  that 
part  of  the  total  activation  energy  needed  to  move  the  impurity  atom  away 
from  its  initial  site  on  plane  1,  (see  Fig.  I),  that  part  needed  to  move  apart 
the  adjacent  lattice  atoms  on  lattice  plane  2  and  thus  permitting  the  impurity 
atom  to  pass  through,  and  E^  the  part  necessary  to  enlarge  the  adjacent 
interstitial  site  on  plane  3  to  accomodate  the  impurity.  Thus  a  very  specific 
distribution  of  the  overall  activation  energy  is  considered  and  it  is  assumed 
that  each  component  is  supplied  at  the  local  temperature  of  the  particular 
plane  involved.  Using  Eq.  (4.12)  we  find  that  the  Wirtz  model  predicts  a 
lattice  bias  W  given  by 


NlTl>  *  ^2  |  “2E2«*3E3  I 


(5.2) 


where  j  an(*  a3  =  exPecte(*»  Ej  does  not  contribute  to  W.  From 

Eqs.  (2.21)  and  (4.21)  the  predicted  heat  of  transport  is  then 

Qm  =  AE  -(2a2E2  *  2a3E3  )  (5*3)‘ 

Wirtz  takes  the  activation  energy  AE  to  be  AE=Ej+E2+Ej  which  yields 


VE3 


(5.4) 


for  the  heat  of  transport.  Equation  (5.4)  predicts  that  the  heat  of  transport 
is  always  less  than  the  activation  energy  AE  and  indeed  obeys  the  stronger 
inequality 


IqJ  ^  ae 


(5.5). 
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This  latter  result  has  been  criticized  since  negative  heats  of  transport, 

(7) 

larger  in  magnitude  than  All,  have  been  observed  experimentally. 

Our  random  walk  treatment,  together  with  the  Einstein  model  of 
independent  oscillators,  has  lead  to  an  equation  for  [Eq.  (4.22)]  which 
is  very  similar  to  Eq.  (5.3)  of  the  Wirtz  model.  The  random  walk  model 
predicts  that  the  heat  of  transport  will  be  less  than  the  activation  energy 
AE,  though  (depending  on  the  values  of  e^,  and  (^negative  heats  of  trans¬ 
port  larger  in  magnitude  than  AE  may  occur.  This  conclusion  is  in  agreement 
with  the  experimental  data  of  Allnatt  and  Chadwick and  of  Oriani.^b) 

The  parameters  E^  and  E^  in  the  Wirtz  model  are  related  to  the  probabilities 
of  achieving  critical  configurations  in  the  host  lattice.  The  cell  model  sug¬ 
gests  that  Ej,  the  part  of  the  activation  enei'gy  related  to  enlarging  the 

3 

adjacent  interstitial  site  on  plane  3,  should  involve  an  of  y  rather  than  1 
since  host  lattice  atoms  must  move  on  the  two  planes  adjacent  to  the  new  inter¬ 
stitial  site.  In  this  case  Eq.  (5.4)  now  reads 


<4  *  V2E3 


(5.6) 


and  the  Wirtz  model  can  predict  a  negative  with  Iq^J  >  AE.  The  Wirtz 
model  also  predicts  that  the  heat  of  transport  is  independent  of  temperature 
(to  the  degree  that  Ej  and  E^  are) ,  though  a  temperature  dependence  is  sometimes 
observed  experimentally.  Both  the  general  random  walk  model,  Eq.  (2.21),  and 
the  cell  model,  Eq.  (4.13),  clearly  permit  a  temperature  dependent  Q^.  When 
the  Einstein  model  approximations  are  valid  and  Eq.  (4.22)  holds,  this  pre¬ 
dicted  temperature  dependence  will  be  small. 

One  can  thus  interpret  the  Wirtz  model  of  thermal  diffusion  as  a  local 
equilibrium  generalization  of  Rice's  isothermal  diffusion  model.  The  specific 
equations  of  the  model,  e.g.  Eq.  (5.1)  for  the  jump  rate  and  Eq.  (5.4)  for 

are  over-simplified  versions  of  a  more  complicated  cell  model  approach,  but 
under  certain  circumstances  the  cell  model  theory  gives  expressions  closely 
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resembling  the  IVirtz  model.  The  physical  picture  behind  the  Wirtz  model 
appears  to  be  a  plausible  and  suggestive  description  of  thermal  diffusion  by 
an  interstitial  mechanism  for  those  systems  where  the  bias  in  jump  direction 
is  primarily  due  to  an  "excluded  volume"  effects. 
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FIGURE  CAPTION 


Figure  1 : 


Schematic  representation  of  interstitial  planes  1  and  3  separated 
by  lattice  plane  2  in  the  one  dimensional  model. 
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